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A viscoplastic theory is developed that reduces analytically to creep theory under 
steady-state conditions. A viscoplastic model is constructed within this theoretical 
framework by defining material functions that have close ties to the physics of 
inelasticity . Asa consequence , this model is easily characterized — only steady-state 
creep data , monotonic stress-strain curves , and saturated stress-strain hysteresis loops 
are required. The model is applied to the copper alloy NARloy Z. 


1 Introduction 

In mankind’s enduring pursuit to go faster and further with 
greater economy and safety in its diverse variety of vehicles 
that travel across land and sea or through the air and space, 
we are taxing our materials to their utmost capabilities. Con- 
sequently, the need for accurate material models to describe 
the various physical properties of a given material is much 
more critical in the design and development of these vehicles 
than it has ever been, and this need can only be expected to 
continue to grow. 

The analysis of metallic response for high temperature ap- 
plications requires mathematical models capable of predicting 
accurately the short-term plastic strains, the long-term creep 
strains, and interactions between them. Viscoplastic models 
attempt to do that. Multiaxial, cyclic and nonisothermal his- 
tories are normal service conditions, not exceptional ones, all 
of which challenge the predictive capabilities of such models. 

Prior to the advent of the computer, viscoplasticity was a 
theory in its infancy; however, over the past two decades sub- 
stantial advancements have been made to the theory. Because 
of viscoplasticity’s innate nature, which leads to systems of 
first-order, ordinary, differential equations that are nonlinear, 
coupled, and mathematically stiff, a unique mathematical 
structure (like that of elasticity) is not to be expected. Never- 
theless, these past two decades have given the community a 
vast wealth of experience with a variety of evolution equa- 
tions — what works, what does not, and in many cases, some 
physical insight as to why. Using this experience base, we have 
set out to develop a viscoplastic model whose predictive ca- 
pabilities are in reasonable agreement with experiments. A 
special emphasis in this development process was that the re- 
sulting model must be characterized easily. The need to follow 
this requirement is vital. Experience has taught us that a mod- 
el’s ease of characterization without calibration via exotic ex- 
periments is often considered by many industrial users of 
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viscoplasticity to be of greater value than the model’s ability 
to predict accurately a material’s behavior (within reason). We 
therefore seek to strike a balance between accuracy and ease 
of characterization, using physics as our guidepost. 

The paper begins with a brief overview of the theories of 
elasticity and creep. This is followed by a definition of vis- 
coplastic flow and the introduction of the required internal 
state variables. The next section demonstrates how a visco- 
plastic theory can be constructed to reduce analytically to creep 
theory under steady-state conditions. This important section 
demonstrates how a bridge between these two theories can be 
built — a concept that is not prevalent in the viscoplastic lit- 
erature. By building this bridge, the model not only has a 
stronger physical base, but it also reduces substantially the 
complexity of material characterization. A succinct description 
of the viscoplastic model is given, and for illustrative purposes, 
the copper alloy NARloy Z is modeled. This material finds 
applications where moderate strength is required under con- 
ditions of very high heat flux, e.g., it is used as the nozzle liner 
material in the main rocket engines of NASA’s space shuttles 
where steep, rapidly applied, thermal gradients cause large 
localized strains. 

2 Elasticity 

The stress, cr,,, is taken to be related to the infinitesimal 
strain, e u , through the constitutive equations of an isotropic 
Hookean material, viz. 

S ij = 2 l i(E ij -efj) where e p kk = 0, (1) 

and 

<?kk = 3x(e** - a ( T— T 0 )b kk )> (2) 

which are characterized by the shear, and bulk, «, elastic 
moduli, and where 

Sy = Ofj - 1/3 o kk b {J and E v = «, y - 1 /3e kk b lJ (3) 

denote the deviatoric stress and strain, respectively. The mean 
coefficient of thermal expansion, ct, acts on the difference 
between the current temperature, T, and some reference tem- 
perature, T 0 . The Kronecker delta, 6 , ; , has the value 1 if /= j t 
otherwise it is 0. Repeated Latin indices are summed from 1 
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Table 1 Elastic constants for NARloy Z (Anonymous, 1986) 


Constants 

Units 

Value 

Of 

K- 1 

16.5 x 10~ 6 

Me 

MPa 

52,000 

Mi 

MPa/K 

-14 

V 

— 

0.34 

M = (i 0 + fx i T, T is in K 


to 3 in the usual manner. Equation (1) characterizes the de- 
viatoric stress response, while Eq. (2) characterizes the hydro- 
static stress response. The plastic strain, e? y , and thermal strain, 
a(T-T 0 )djj 9 are, in essence, eigenstrains that represent de- 
viations from deviatoric and hydrostatic elastic behaviors, re- 
spectively. 

Young’s modulus, E, and Poisson’s ratio, v , are the two 
elastic constants that are usually determined via experiment. 
The expressions, 


E 

2(1 + v) 


and 


E 

30-2*)’ 


(4) 


define their interdependence with the elastic moduli of Eqs. 
(1) and (2). Only two elastic moduli are independent for elast- 
ically isotropic materials. Values for the elastic constants of 
NARloy Z (typical composition: Cu-3%Ag-0.5%Zr) are given 
in Table 1 . 


3 Creep 

The evolution of plastic strain which describes the classical 
theory of creep (Odqvist, 1974) is given by 

ef y l H =l/2«e p ll„^, (5) 

with the subscript “ss” implying steady-state, and where 
le p ll denotes the magnitude of plastic strain-rate, and IIS II 
denotes the magnitude of deviatoric stress. The subscript sig- 
nifying steady state is not attached to Sjj because stress is a 
controllable external variable, whereas creep rate is a response 
variable. This equation states that an increment in creep strain 
accumulates in the current direction of the deviatoric stress. 
A dot is placed over a variable to signify its time rate-of-change. 

The norms, or magnitudes, pertaining to the deviatoric ten- 
sors of this paper are defined by 

III I =ypn~hj and (6) 


where I tJ is any deviatoric “strain-like” tensor, and J tj is any 
deviatoric “stress-like” tensor. These are the norms of von 
Mises (1913), where the coefficients under the radical signs 
scale the theory for shear. 

In the theory of creep, I e p ll^ is described by a kinetic equa- 
tion, i.e., an equation of state. Zener and Hollomon 0 944) 
determined that such a kinetic equation can, to a good ap- 
proximation, be decomposed into a product of two functions; 
in particular, at steady state 


U p l„ = fi[T\Z„ 



(7) 


where 0 >0 is a thermal function, Z> 0 is the Zener parameter, 
and C> 0 is a strength parameter that normalizes the stress. 
The Zener parameter is a temperature normalized measure of 
the plastic strain-rate. Square brackets, [•], are used through- 
out this paper to denote “function of,” and are therefore kept 
logically separate from parentheses, (•), which are used for 
mathematical groupings. 

In the physical description of the thermal function, #, there 
is a parameter called the activation energy, Q y which — for creep 
at low stresses and elevated temperatures — is associated with 
self-diffusion where the rate-controlling mechanism for de- 
formation is dislocation climb (Sherby and Weertman, 1979). 


At higher stresses and/or more moderate temperatures, the 
rate-controlling mechanism changes from diffusion-controlled 
dislocation climb to obstacle-controlled dislocation glide (Kocks 
et al., 1975). Along with this change in the deformation mech- 
anism, there occurs a change in the activation energy (Sherby 
and Burke, 1968). Miller (1976) approximates the observed 
temperature dependence of the activation energy for steady- 
state flow with a linear function for temperatures below some 
threshold temperature, T ly while for temperatures above this 
threshold the activation energy is kept constant, in accordance 
with the experimental observations of Dorn (1954). Because it 
is the free energy (not the activation energy) that drives the 
kinetics of plastic deformation (Kocks et al., 1975), Miller 
integrated his linear function for the activation energy and 
obtained the following Arrhenius-like expression for the ther- 
mal function, 



when 7 f < T<T m 

( 8 ) 

when 0 <T<T, 


where k is the universal gas constant (8.314 J/mole*K). The 
applicability of this relationship is discussed elsewhere (Freed 
et al., 1992). The transition temperature, T { , between these 
two domains in activation energy is not unique; it is known 
to depend on the strain-rates used to make the measurements 
for activation energy. An increase in strain-rate increases the 
transition temperature (Sherby and Burke, 1968). For the vast 
majority of engineering applications, a transition temperature 
of r r = 1/2 T m seems appropriate for f.c.c. metals, and is used 
in our characterization of NARloy Z. 

When the mechanism for deformation changes from dif- 
fusion-controlled dislocation climb to obstacle-controlled dis- 
location glide, the creep response changes from power-law to 
exponential behavior (Ashby, 1972). Following the approach 
of Miller (1976), we adopt Garofalo’s (1963) empirical expres- 
sion for the steady-state Zener parameter, i.e., 



where ^>0, C> 0 and rt>0 are the material constants. For 
stress states below power-law breakdown, i.e., when IIS I! <C, 
the steady-state Zener parameter of Garofalo reduces to the 
power-law relationship 

/iisiiy 

Zss=aI— j , do) 


thereby designating dislocation climb as the rate-controlling 
mechanism. (Note: A , C, and n are independent in Eq. (9) but 
not in Eq. (10).) Similarly, when the stress exceeds powerlaw 
breakdown, i.e., when IISil>C, Garofalo’s Zener parameter 
reduces to the exponential relationship 


Z SS = A ' exp 


IISI! 

C' ’ 


( 11 ) 


where A ' =,4/2" and C' -C/n y thereby designating disloca- 
tion glide as the rate-controlling mechanism. The ability of 
Eqs. (7)-(9) to correlate the stationary creep-rate data of NAR- 
loy Z is demonstrated in Fig. 1 . The material constants obtained 
from this correlation are given in Table 2. Because none of 
these data lie within the power-law domain, the exponential 
creep equation, Eq. (11), was used to determine values for A ' 
and C' leading to the straight line fit shown in the log/linear 
plot of Fig. 1(a), where A ' = 5 x 10 1 ' s" 1 and C' - 3.5 MPa 
for the predefined values of Q = 450,000 J/mole (Lewis, 1970) 
and 7> 400°C (assumed). 1 Taking a = 4 (assumed), the values 


1 This seems to be an excessively large value for Q , but it is the only experi- 
mentally determined value currently available to us. 
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Fig. 1(a) 



Fig. 1(6) 

Fig. 1 Steady state creep behavior of NARloy Z. Data are from Lewis 
(1970) with the reported material composition: Cu-2.89%Ag-0.22%Zr. 


for A ' and C' have been converted to those of A and C that 
are reported in Table 2. The result is the curved line presented 
in the log/log plot of Fig. 1(6). We note that the value of C 
for NARloy Z, i.e., 14 MPa, obtained with this choice for n, 
i.e., 4, is in agreement with the value of C for Cu, i.e., 13 
MPa, reported in Freed and Walker (1993a). 

This continuum representation for creep is well established. 
Our viscoplastic model reduces analytically to this creep model 
under steady-state conditions. Hence, the material constants 
that characterize this creep model also appear in our visco- 
plastic model, which simplifies substantially its characteriza- 
tion process. 

4 Viscoplastic Flow 

A general mathematical structure for viscoplasticity (Freed 
et al., 1991) may admit up to three kinds of internal state 
variables; they are: (i) the (scalar-valued) drag strength, D> 0; 
(ii) the (scalar-valued) yield stress, Y> 0; and (iii) the (devia- 
toric tensor-valued) back stress, B tJ . The drag strength and 
yield stress account for isotropic hardening effects, while the 
back stress accounts for kinematic (flow-induced anisotropic) 
hardening effects. 


Table 2 Steady-state creep constants for NARloy Z 


Constant 

Units 

Value 

A 

s" 1 

8 x 10 18 

C 

MPa 

14 

rt 

— 

4 

Q 

J/moi. 

450,000 

T m 

K 

- 1350 

T,= l/2T m 


Prager’s (1949) constitutive relation is used to describe the 
evolution of viscoplastic flow-, i.e., 




1/211 t p 


Sg Bjj 

II S — B II ’ 


( 12 ) 


This particular choice for the flow law implies that a nested 
set of flow surfaces exists; they are surfaces of constant plastic 
strain-rate when evaluated under isothermal conditions. This 
constitutes a set of ellipsoids in deviatoric stress space that are 
centered on the back stress. 

The kinetics of viscoplasticity are taken to be described by 
a Zener and Hollomon (1944) type decomposition of state 
(Freed et ah, 1991), viz. 


\\t p \\ = d[T)Z 


^IIS-B> - Y \ 


* 0 , 


(13) 


where the Macauley bracket, <(IIS- BH - Y)/D), has either a 
value of 0 whenever IIS-BH < Y (defining the elastic domain), 
or a value of ( II S - B II - Y ) /D whenever H S - B II > Y (de f i ning 
the viscoplastic domain), with IIS - B II = Y establishing the yield 
surface. Many viscoplastic models have no distinct yield sur- 
face, i.e., they set Y to 0. The distinguishing feature between 
viscoplasticity (a rate-dependent theory) and plasticity (a rate- 
independent theory) is that viscoplasticity admits states both 
inside and outside of the yield surface (governed by a kinetic 
equation of state); whereas, plasticity admits only states that 
are inside and on the yield surface (governed by a consistency 
condition), but not outside of it. As a consequence, the plastic 
strain-rate is continuous as one moves from the elastic domain 
across the yield surface and into the inelastic domain of vis- 
coplastic response; whereas, the accumulation of plastic strain 
is discontinuous as one moves from the elastic domain onto 
the yield surface in plasticity. The elastic domain of many 
viscoplastic models is shrunk to a point, as they do not admit 
a yield surface. 

The internal state variables — B ijy D and Y — are described 
by evolution equations that are functions of state. The back 
stress evolves rapidly when compared with the rates of evo- 
lution for the drag strength and yield stress, which is a source 
of mathematical stiffness in the governing equations of vis- 
coplasticity. The evolution of the back stress accounts for the 
change in material stiffness that is observed during the tran- 
sition from elastic to plastic behavior, while the evolutions of 
the drag strength and yield stress account for the more gradual 
work hardening processes that are caused by the overall ac- 
cumulation of plastic deformation. The internal variables are 
considered to evolve phenomenologically through competitive 
processes associated with strain hardening, strain-induced dy- 
namic recovery, and time-induced thermal recovery. Their spe- 
cific functional forms are presented later in Section 6, whose 
derivations are given in the conference proceedings’ version 
of this paper, i.e., Freed and Walker (1993b). 

5 Creep <=* Viscoplasticity 

In the process of going from creep theory to viscoplasticity, 
one must remove the steady-state constraint that is present in 
creep, and thereby extend the domain of admissible states to 
include transient behavior. In other words, viscoplasticity is 
capable of modeling both primary and secondary creep be- 
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havior. The modeling of transient behavior is done through 
the introduction of internal state variables. Although the pur- 
pose of viscoplasticity is to model rate-dependent transient 
behavior, it is not unreasonable to also require that it reduces 
to creep theory under steady-state conditions. An important 
objective in our development of a viscoplastic theory is that 
it reduces analytically to creep theory when at steady state. 
Not only is this a realistic requirement, but it also strengthens 
the physics of the theory, and it simplifies greatly the process 
of model characterization — about half of our viscoplastic ma- 
terial constants come from correlating stationary creep-rate 
data alone. 

In order for a viscoplastic theory to reduce analytically to 
creep theory when at steady-state (i.e., when 8 = 0, D = 0, and 
7=0 for I e p ll >0) two conditions must be satisfied. First, the 
back stress must be coaxial with the stress at steady state so 
that the directions of plastic strain-rate defined by Eqs. (5) 
and (12) are also coaxial at steady state. And second, it is 
necessary that the kinetics of viscoplasticity, Eq. (13), reduce 
analytically to the kinetics of creep, Eq. (7), under steady-state 
conditions. The evolution law for back stress given in Eq. (28) 
satisfies this first constraint. To satisfy the second constraint, 
one must first hypothesize a relationship between the steady- 
state and transient Zener parameters, and then hypothesize 
another one between the internal and external variables, when 
at steady state (Freed and Walker, 1990). We therefore suppose 
that 




IIS -B II 
D 



(14) 


in support of Eq. (13). This relationship implies that the tran- 
sient Zener parameter, Z, has the same functional form as the 
steady-state Zener parameter, Z ss > but with a different argu- 
ment; in particular, and in accordance with Eq. (9), we take 


Z = A sinh* 


IIS-BII 



(15) 


which is similar in form to the kinetics of Miller's (1976) vis- 
coplastic model, but with a yield stress and without a power 
acting on the Macauley bracket. 

Furthermore, we shall suppose that 

IIBIU=/UIISII]IISII, A, = A> + *HSII 

and y’„«(l-/)t Jf [lSHISI, (16) 

in support of experimental evidence, where t 5S >0 and 8 >0 are 
the steady-state fractions of applied stress that are associated 
with the internal stress (i.e., the back and yield stresses) and 
the drag strength, respectively, such that \/2<i ss < 1. The pa- 
rameter / partitions the internal stress between isotropic and 
kinematic contributions, such that 0</< 1. The fact the drag 
strength is taken to be proportional to the saturation stress is 
a consequence of the fact that the drag strength represents the 
material’s innate strength to resist plastic flow, i.e., D is a 
strength parameter — not a stress parameter. We take the in- 
ternal stress to be a nonlinear function of the applied stress at 
saturation because that is what the experimental data of Argon 
and Takeuchi (1981) and Cadek (1987) suggest. A similar hy- 
pothesis to that of Eq. (16) is given in Freed and Walker (1993a) 
for the case where the internal stress is composed of two back 
stresses with no yield stress. 

Because the applied stress and the back stress must be coaxial 
at steady state, as discussed above, it follows that 

IIS-BI*=IS1- IBI«. (17) 

Therefore, upon equating the arguments of the Zener param- 
eters in Eqs. (7) and (14), while utilizing Eqs. (16) and (17), 
one obtains the result 


E0.75 h 


"0.50 


of 0.25 


0.00 



D n 


0.00 


0.25 


0.50 

ISI/ISL 


0.75 


1.00 


Fig. 2 Schematic of internal stress (back plus yield stresses) at steady 
state versus applied stress 


c-Do-smw 


C-£> 0 + V (C-Z?q) 2 - 46C( HBHjj-h P c ) 
2C 


(18) 


If one uses Eq. (16) and writes HBII W + y^=t„[lSH]iSII, then 
from Eq. (18) one determines that 11811^ + = (C-Z) 0 

- 8 IIS II) IIS II /C. Because 3(181^ + T^/dllSIl = 0 establishes the 
maximum state of internal stress, one is lead to the result 

»Sll max = ^p. (19) 


Substituting this relationship back into Eqs. (16) and (18) gives 
additional upper bounds for: the back stress, 



(20) 

the drag strength, 


Anax= 1/2(C+Z) 0 )> 

(21) 

and the yield stress, 


, =(1 -n 

max J) 4 8C ‘ 

(22) 


Similar bounds are given in Freed and Walker (1993a) for the 
case where the internal stress is composed of two back stresses 
and no yield stress. It is a remarkable fact that one can bound 
the stress and internal state variables without specifying any- 
thing about how these internal state variables evolve . 

Restricting i S5 to be real valued, and considering ^ to be 
associated with the maximum attainable magnitude of internal 
stress, one finds on approaching the limit of zero stress that 
the ratio of internal stress to applied stress at steady state is 
at its maximum, i.e., 


, _ C-A> 

lim Ljj — l m ax — 
isi-o C 


(23) 


which is in reasonable agreement with Argon and Takeuchi’ s 
(1981) and Cadek’s (1987) experimental observations. Ap- 
proaching the limit of maximum stress, this ratio attains its 
minimum, i.e., 


lim t aStrain = %^=l/2, (24) 

IISI — SSImax /c 


which is in reasonable agreement with Lowe and Miller’s ( 1 983) 
and Argon and Bhattacharya’s (1987) experimental observa- 
tions. A schematic of the steady-state internal stress versus the 
applied stress — as predicted by Eqs. (16) and (18) with typical 
values of Do -C / 100 and /= 0.6 — is presented in Fig. 2. The 
trends depicted therein are in qualitative agreement with the 
experimental results referenced above. 
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Table 3 Additional viscoplastic constants for NARloyZ 


Constants 

Units 

Value 

6 



0.035 

f 



0.65 

b 

MPa 

11 

t 

— 

6.4 

D 0 = C/100 


To be physically meaningful, II B it >0, Z)>0, and Y> 0. Fur- 
thermore, their steady-state values ought to increase mono- 
tonically with increasing stress (Freed and Walker, 1990). This 
is verified easily for our hypothesis, Eqs. ((14), (16), and (17)), 
as long as 0< IIBII< IIBII max , D 0 <D<D max and 0< Y — ^max- 



6 The Model 

A succinct description of our viscoplastic model is given 
below. The stress is acquired through the constitutive equations 

S IJ = 2^(E tJ - and <J kk = 3tc(e kk - a{T- T 0 )b kk ) . (25) 

The flow equation and kinetics that describe plastic straining 
are given by 

kfj= 1/211 * P| isrii and un=dz ' (26) 


respectively, with the von Mises norm of effective stress being 
defined by 

IIS- Bll = \ 1 /2(Sij- J~j ) ( S„ (27) 
The evolutions of back stress and drag strength are given by 

a, - 2 «(<r, 

and £> = /t(ll^ll-Alle p ll-t?r), (28) 


respectively, such that Z? 0 ^7)<Z) max , while the yield stress is 
related through the state function 


Y=(\-f) 


( D-D 0 )(C-D ) 
bC 


(29) 


which is not an evolution equation. Associated with these re- 
lationships are the material functions: 




-Q 


exp 

' 

kT 


| 

iO 

i 

' exp 

I 

1 

1C 


Z — / 


In 


A 

+ oi 

T 

/J 


when T t <T<T m 
when 0< T, y 

IIS -Bll - Y\ 

D 


(30) 


(31) 


H=(0A+0.9£)ti and 


{D- D 0 ) (C-D) 
} bC 


(32) 


with 


A = f£ 


and 


r-A sinh" 


D-Do 

bC 


(33) 




1 ( ^ ! J / Sjj _ Bjj \ 

\8\HS-BII L J y IIS — Bll L )' 


(34) 


Restricting the drag strength to be bound by the interval 
D Q <D<D m ^ restricts automatically the remaining variables: 
0< IlSil < IISIl max , 0< II Bll < IIB!l max and 0< Y< Y max . The de- 
velopment of the evolution equations for the back stress and 
drag strength, along with the derivations of their associated 
material functions, are given in the conference proceedings’ 
version of this paper (Freed and Walker, 1993b). Also found 
therein is a detailed discussion of how one goes about char- 


Fig. 3 Tensile stress-strain curves for NARloy Z. Data are from the 
Anonymous (1986). Experiment ( o) versus theory (—). it =0.001 s'* 1 (as- 
sumed), Dl,. o = 2.02 MPa (for T= 21 and 538°C) and Dl„ 0 =2.64 MPa 
(for T=315 and 427°C). 



£ 


Fig. 4 Saturated hysteresis loops for NARloy Z. Data are from Conway 
et al. (1975). Experiment (o) versus theory (— ). « =0.004 s 1 and 
T=538°C. 


acterizing this particular model. Values for the remaining tran- 
sient constants of our viscoplastic model are given in Table 3 
for NARloy Z. 

For unaxial loading histories in tension and compression, 
the above governing equations hold with the following alter- 
ations: 


o = E(e-e p -a(T- T 0 )), 


e p = sgn[a-/3] 


II € ^ II 

A 


and 


l = 3//(e p -^- || 6 P|1 ) 


(35) 


given that a = an = 3/2S n , 0- 0 U = 3/2 B ih e-e n and z p = e p \. 
As for the material functions, the above equations apply with 
the following alterations: 


Z — A sinh n 


ig-gi zAh Y 

Ad /_ 


and ?=l/2 


sgn[a-/3]l 


Al 


(36) 


The ability of the model to correlate (not predict) monotonic 
and cyclic material behavior is presented in Figs. 3 and 4, 
respectively, for NARloy Z. Data for this material are sparse, 
thereby not permitting a more detailed assessment of the mod- 
el’s predictive capability. In part, this demonstrates a design 
objective in our development of this model — the capability to 
characterize the model from a sparse data set. 

This model is not perfect, and certainly not ideal, but hope- 
fully it represents another step in that direction. It is a sim- 
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plified continuum description of complex microscopic 
phenomena; nevertheless, its development has been guided by 
the physics of these phenomena. There are several known de- 
ficiencies associated with this model. They are; the predicted, 
transient, rate dependence, which is extrapolated from steady- 
state dependence, does not always match experimentally ob- 
served rate dependence, for example, the rate dependence ex- 
hibited during a stress relaxation experiment (Freed and Walker, 
1993a); predicted transient behavior in the region of transition 
between the domains of power-law and exponential behaviors, 
which is also taken from steady-state behavior, does not always 
agree with experimental observations (Loh, 1993); and the well- 
known fact of excessive, predicted, ratchetting behavior, which 
is a consequence of our using the Armstrong and Frederick 
(1966) evolution equation for back stress (Freed and Walker, 
1993c). 


7 Closure 

By designing the development of our viscoplastic model in 
such a manner that it reduces analytically to a creep model 
under steady-state conditions, we have incorporated essential 
physics into our model, and we have also simplified greatly 
the process that one must go through in order to completely 
characterize a material with this model. In this sense, we have 
developed a viscoplastic model with an eye towards its char- 
acterization. This has particular merit because parameter es- 
timation of a viscoplastic model is, in general, a very complex 
process that all too often prohibits its use in applications. A 
model’s relative ease of characterization without the need for 
exotic experiments is often considered by many industrial users 
of viscoplasticity to be of greater value than the model’s ability 
to predict accurately a material’s behavior (within reason). Our 
model was developed with this fact in mind, where we have 
sought to strike a balance between accuracy and ease of char- 
acterization using physics as our guidepost. 
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